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Abstract 

We perform a detailed study of a variety of LHC signals in supersymmetric models 
where lepton number is promoted to an (approximate) symmetry. Such a sym- 

metry has interesting implications for naturalness, as well as flavor- and CP-violation, 
among others. Interestingly, it makes large sneutrino vacuum expectation values phe- 
nomenologically viable, so that a slepton doublet can play the role of the down-type 
Higgs. As a result, (some of) the leptons and neutrinos are incorporated into the 
chargino and neutralino sectors. This leads to characteristic decay patterns that can 
be experimentally tested at the LHC. The corresponding collider phenomenology is 
largely determined by the new approximately conserved quantum number, which is 
itself closely tied to the presence of "leptonic R-parity violation" . We find rather loose 
bounds on the first and second generation squarks, arising from a combination of sup- 
pressed production rates together with relatively small signal efficiencies of the current 
searches. Naturalness would indicate that such a framework should be discovered in 
the near future, perhaps through spectacular signals exhibiting the lepto-quark nature 
of the third generation squarks. The presence of fully visible decays, in addition to de- 
cay chains involving large missing energy (in the form of neutrinos) could give handles 
to access the details of the spectrum of new particles, if excesses over SM background 
were to be observed. The scale of neutrino masses is intimately tied to the source of 
[/(I) /J breaking, thus opening a window into the i?- breaking sector through neutrino 
physics. Further theoretical aspects of the model have been presented in the companion 
paper [ij. 
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1 Introduction 



The recent discovery at the LHC of a Higgs-hke signal at ~ 125 GeV has put the general issue 
of electroweak symmetry breaking under a renewed perspective. In addition, the absence of 
other new physics signals is rapidly constraining a number of theoretically well-motivated 
scenarios. One of these concerns supersymmetry, which in its minimal version is being 
tested already above the TeV scale. In view of this, it is pertinent to consider alternate 
realizations that could allow our prejudices regarding e.g. naturalness to be consistent with 
the current experimental landscape, within a supersymmetric framework. At the same time, 
such scenarios might motivate studies for non-standard new physics signals. 

One such non-standard realization of supersymmetry involves the possible existence of 
an approximately conserved i?-symmetry at the electroweak scale [2 -13 . It is known that 



one of the characteristics of such scenarios, namely the Dirac character of the gauginos 



(in particular, gluinos), can significantly soften the current exclusion bounds 14,15 . At 
the same time, an approximate /^-symmetry which extends to the matter sector, could end 
up playing a role akin to the GIM mechanism in the SM, thereby allowing to understand 



the observed flavor properties of the light (SM) particles. As advocated in Ref. 16 , a 
particularly interesting possibility is that the i?-symmetry be an extension of lepton number 
(see also [l7]). In a companion paper ^J, we classify the phenomenologically viable R- 
symmetric models, and present a number of theoretical and phenomenological aspects of 
the case in which R-symmetry is tied to the lepton number. Such a realization involves 
the "R-parity violating (RPV) superpotential operators" , W D XLLE'^ + X'LQD'^ where, 
unlike in standard RPV scenarios, there is a well-motivated structure for the new A and A' 
couplings, some of them being related to (essentially) known Yukawa couplings. Although, 
at first glance, one might think that such a setup, possibly with a preponderance of leptonic 
signals should be rather constrained, we shall establish here that this is not the case. In fact, 
the scenario is easily consistent with most of the superparticles lying below the TeV scale. 
Only the Dirac gauginos are expected to be somewhat above the TeV scale, which may be 
completely consistent with naturalness considerations. As we will see, the light spectrum is 
particularly simple: there is no LR mixing in the scalar sector, and there is only one light 
(Higgsino-like) neutralino/chargino pair. At the same time, it turns out that the (electron) 
sneutrino vev can be sizable, since it is not constrained by neutrino masses (in contrast to 
that in standard RPV models). This is because the Lagrangian (approximately) respects 
lepton number, which is here an R symmetry, and the sneutrinos do not carry lepton number. 
Such a sizable vev leads to a mixing of the neutralinos/charginos above with the neutrino 
and charged lepton sectors (z/g and e~ to be precise), which results in novel signatures and 
a rather rich phenomenology. Although the flavor physics can in principle also be very rich, 
we will not consider this angle here. 

We give a self-contained summary of all the important physics aspects that are relevant 
to the collider phenomenology in Section |2} This will also serve to motivate the specific 
spectrum that will be used as a basis for our study. In Section [3| we put together all the 
relevant decay widths, as a preliminary step for exploring the collider phenomenology. In 
Section |4] we discuss the current constraints pertaining to the first and second generation 
squarks, concluding that they can be as light as 500 — 700 GeV. We turn our attention to the 
third generation phenomenology in Section |5| where we show that naturalness considerations 
would indicate that interesting signals could be imminent, if this scenario is relevant to the 
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weak scale. In Section |6| we summarize the most important points, and discuss a number 
of experimental handles that could allow to establish the presence of a leptonic R symmetry 
at the TeV scale. 



2 U{1)r Lepton Number: General Properties 

Our basic assumption is that the Lagrangian at the TeV scale is approximately U{1)r sym- 
metric, with the scale of U{1)r symmetry breaking being negligible for the purpose of the 
phenomenology at colliders. Therefore, we will concentrate on the exact /^-symmetric limit, 
which means that the patterns of production and decays are controlled by a new (approx- 
imately) conserved quantum number. We will focus on the novel case in which the R- 
symmetry is an extension of the SM lepton number. Note that this means that the extension 
of lepton number to the new (supersymmetric) sector is non-standard. 



2.1 The Fermionic Sector 

As in the MSSM, the new fermionic sector is naturally divided into strongly interacting 
fermions (gluinos), weakly interacting but electrically charged fermions (charginos) and 
weakly interacting neutral fermions (neutralinos) . However, in our framework there are 
important new ingredients, and it is worth summarizing the physical field content. This will 
also give us the opportunity to introduce useful notation. 



2.1.1 Gluinos 

One of the important characteristics of the setup under study is the Dirac nature of gauginos. 
In the case of the gluon superpartners, this means that there exists a fermionic colored octet 
(arising from a chiral superfield) that marries the fermionic components of the SU{'i)c vector 
superfield through a Dirac mass term: M-f^^o"" + h.c., where a is a color index in the adjoint 
representation of SU{?))c and a is a Lorentz index (in 2-component notation). Whenever 
necessary, we will refer to o as the octetino components, and to g as the gluino components. 
For the most part, we will focus directly on the 4-component fermions = {g'^, o"""} and we 
will refer to them as (Dirac) gluinos, since they play a role analogous to Majorana gluinos in 
the context of the MSSM. However, here the Majorana masses are negligible (we effectively 
set them to zero) and, as a result, the Dirac gluinos carry an approximately conserved (R) 
charge. In particular, R{g) = —R{o) = 1, so that R{G) = 1. -R-charge (approximate) 
conservation plays an important role in the collider phenomenology. 

The Dirac gluino pair-production cross-section is about twice as large as the Majorana 
gluino one, due to the larger number of degrees of freedom. Assuming heavy squarks, and 



within a variety of simplified model scenarios, both ATLAS 18 20 and CMS 21 24 have 
set limits on Majorana gluinos in the 0.9 — 1 TeV range. As computed with Prospino2 pS] 
in this limit of decoupled squarks, the NLO Majorana gluino pair-production cross-section is 
"^Majorana ("^9 = 1 TcV) ~ 8 fb at the 7 TcV LHC run. Although, for the same mass, the Dirac 
gluino production cross-section is significantly larger, it also falls very fast with the gluino 
mass so that the above limits, when interpreted in the Dirac gluino context, do not change 
qualitatively. Indeed, assuming a similar K-factor in the Dirac gluino case, we find a NLO 
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pair-production cross section of 

(^dLc(^3^ = 1-08 TeV) ^ 8 fb. Nevertheless, from a theo- 
retical point of view the restrictions on Dirac gluinos coming from naturalness considerations 
are different from those on Majorana gluinos, and allow them to be significantly heavier. We 
will take Mg = = 2 TeV to emphasize this aspect. This is sufficiently heavy that direct 
gluino pair-production will play a negligible role in this studyj^ At the same time, such 
gluinos can still affect the pair-production of squarks through gluino t-channel diagrams, as 
discussed later (for the gluinos to be effectively decoupled, as assumed in e.g. [15j, they must 
be heavier than about 5 TeV). 

2.1.2 Charginos 

We move next to the chargino sector. This includes the charged fermionic SU (2)^ superpart- 
ners (winos) and the charged tripletino components, Tj" and , of a fermionic adjoint 
of SU{2)l (arising from a triplet chiral superfield). It also includes the charged components 
of the Higgsinos, and r^. The use of the notation instead of h'^ indicates that, unlike 
in the MSSM, the neutral "Higgs" component i?^ does not acquire a vev. Rather, in our 
setup, the role of the down-type Higgs is played by the electron sneutrino Vf. (we will denote 
its vev by Vg). As a result, the LH electron mixes with the above charged fermions, and 
becomes part of the chargino sector (as does the RH electron field e^). Besides the gauge 
interactions, an important role is played by the superpotential operator W D X^HuTRa, 
where T is the SU{2)l triplet superfield [I]. 

The pattern of mixings among these fermions is dictated by the conservation of the 
electric as well as the _R-charges: i?(w^) = -R(e^) = R{r~^) = +1 and R(T^) = R(T^) = 
R{e2) = R{h^) = —1. In 2-component notation, we then have that the physical charginos 
have the composition 
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where i = 1,2. The notation here emphasizes the conserved electric and i?-charges, by 
indicating them as superindices, e.g. Xi"" denoting the two charginos with Q = +1 and 
R = —1. The U^, are 2x2 unitary matrices that diagonalize the corresponding chargino 
mass matrices. The superindex denotes the product RxQ, while the subindices in the matrix 
elements should have an obvious interpretation. We refer the reader to Ref. for further 
details. In this work we will not consider the possibility of CP violation, and therefore all the 
matrix elements will be taken to be real. The above states are naturally arranged into four 
4-component Dirac fields, X^~^ = {xt~^^Xi ) and = {xt~,x7^)y f*^^ i = ^y'^y whose 
charge conjugates will be denoted by and X^'^. In this notation, e = X^ corresponds 
to the physical electron (Dirac) field. 

^However, at 14 TeV, with o'D^;.a^j.(Af-P ^ 2 TeV) w 3 fb, direct gluino pair-production may become 
interesting. The K-factor (« 2.6) is taken from the Majorana case, as given by Prospino2. This production 
cross-section is dominated by gluon fusion, and is therefore relatively insensitive to the precise squark masses. 
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As explained in the companion paper [Tj, precision measurements of the electron prop- 
erties place bounds on the allowed admixtures and f/^, that result in a lower bound on 
the Dirac masses, written as M^{w~^T^ + w~T^) + h.c. This lower bound can be as low as 
300 GeV for an appropriate choice of the sneutrino vev. However, a sizably interesting range 
for the sneutrino vev requires that be above about 1 TeV. For definiteness, we take in 
this work = 1.5 TeV, which implies that 10 GeV ^ Ve < 60 GeV. Thus, the heaviest 
charginos are the ^ {w^,T^) and ^ (T^,w~) Dirac fields, which we will simply 
call "winos". The lightest chargino is the electron, e ~ (^2^,6^), with non-SM admixtures 
below the 10~^ level. The remaining state is expected to be almost pure hu-Vd, with a mass 
set by the /i-termj^ Naturalness considerations suggest that this parameter should be around 
the EW scale, and we will take /i = 200 — 300 GeV. However, it is important to note that the 
gaugino component of this Higgsino-like state, U{^, although small, should not be neglected. 
This is the case when considering the Xi~ couplings to the first two generations, which cou- 
ple to the Higgsino content only through suppressed Yukawa interactions. In the left panel 
of Fig. [2| we exhibit the mixing angles of the two lightest chargino states as a function of the 
sneutrino vev, Vg, for = 1.5 TeV, fi = 200 GeV, = and = 1. The l^-type matrix 
elements are shown as solid lines, while the ?7-type matrix elements are shown as dashed 
lines (sometimes they overlap). In the right panel we show the chargino composition as a 
function of A^ for Vg = 10 GeV. This illustrates that there can be accidental cancellations, as 
seen for the w~ component of X^~^ at small values of A^. For the most part, we will choose 
parameters that avoid such special points, in order to focus on the "typical" cases. It is 
also important to note that the quantum numbers of these two lightest chargino states (the 
lightest of which is the physical electron) are different. This has important consequences for 
the collider phenomenology, as we will see. 

2.1.3 Neutralinos 

The description of the neutralino sector bears some similarities to the chargino case discussed 
above. In particular, and unlike in the MSSM, it is natural to work in a Dirac basis. The 
gauge eigenstates are the hypercharge superpartner (bino) b, the neutral wino w, a SM 
singlet, s, the neutral tripletino T°, the neutral Higgsinos, and f° and, finally, the electron- 
neutrino z/g (which mixes with the remaining neutralinos when the electron sneutrino gets 
a vev). If there were a right-handed neutrino it would also be naturally incorporated into 
the neutralino sector. In principle, due to the neutrino mixing angles (from the PMNS 
mixing matrix) the other neutrinos also enter in a non-trivial way. However, for the LHC 
phenomenology these mixings can be neglected, which we shall do for simplicity in the 
following. Besides the gauge interactions and the A^ superpotential coupling introduced in 
the previous subsection, there is a second superpotential interaction, W D X^SHuRd, where 
S is the SM singlet superfield, that can sometimes be relevant [l]. 

■^In the companion paper [l] we have denoted this /z-term as /i„ to emphasize that it is different from 
the "standard" /z-term: the former is the coefficient of the H^Rd superpotential operator, where Rd does 
not get a vev and, therefore, does not contribute to fermion masses, while the role of the latter in the 
present scenario is played by /i'iJ„_Le, with being the electron doublet whose sneutrino component gets 
a non-vanishing vev. While the first one is allowed by the U{1)r symmetry, the second one is suppressed. 
However, for notational simplicity, in this paper we will denote the U{l)ii preserving term simply by /x, since 
the "standard", f/(l)i? violating one, will not enter in our discussion. 
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Figure 1: Composition of the two lightest chargino states as a function of the sneutrino vev (left panel) 
and as a function of (right panel). We fix — 1.5 TeV, ^ — 200 GeV and Af = 0. In the left panel we 
take AJ = 1 and in the right panel we take Vg = 10 GeV. We plot the absolute magnitude of the rotation 
matrix elements V^^ (solid lines) and (dashed lines). Not plotted are y+. = and Vi+ 1. X++ is the 
physical (charge conjugated) electron, and is the lightest BSM chargino state (which is Higgsino-like) . 

For reference, we also show in the upper horizontal scale the values of tan/? = v^Jv^. 



In two- component notation, we have neutralino states of definite U{1)r charge 

xr = Uf,s + U^,f' + UZhl + UZ^e, (2) 
where V^^ and are the unitary matrices that diagonahze the neutrahno mass matrix 



(full details are given in Ref. [l 



). These states form Dirac fermions X^^ = (Xi'"''; Xi"); foi' 
i = 1,2,3, where, as explained in the previous subsection, the superindices indicate the 
electric and i?-charges. In addition, there remains a massless Weyl neutralino: 

Xt = Ug~s + U^,f' + Ulhl + UZv^, (3) 

which corresponds to the physical electron- neutrino. With some abuse of notation we will 
refer to X4~ as 'Vg" in subsequent sections, where it will always denote the above mass 
eigenstate and should cause no confusion with the original gauge eigenstate. Similarly, we 
will refer to as the "lightest neutralino" , with the understanding that strictly speaking 
it is the second lightest. Nevertheless, we find it more intuitive to reserve the nomencla- 
ture "neutralino" for the states not yet discovered. The heavier neutralinos are labeled 
accordingly. 

Given that both the gluino and wino states are taken to be above a TeV, we shall also take 
the Dirac bino mass somewhat large, specifically = 1 TeV. This is mostly a simplifying 
assumption, for instance closing squark decay channels into the "second" lightest neutralino 
(which is bino-like). Thus, the lightest (non SM-like) neutralino is Higgsino-like, and is fairly 
degenerate with the lightest (non SM-like) chargino, . 

In Fig. |2| we show the composition of the physical neutrino (x4~) and of the Higgsino- 
like neutralino state {X^'^). Note that, as a result of i?-charge conservation, the neutrino 
state has no wino/bino components. In addition, its (up-type) Higgsino component is rather 
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Figure 2: Left panel: x"" (neutrino) and right panel: (Higgsino-like) composition for — 1 TeV, 

M2 = 1.5 TeV, /X = 200 GeV, A.f = and = 1, as a function of the sneutrino vev. We plot the absolute 
magnitude of the rotation matrix elements and U^. 

suppressed. As a result, the usual gauge or Yukawa induced interactions are very small. 
Instead, the dominant couplings of X4~ to other states will be those inherited from the 
neutrino content itself. The associated missing energy signals will then have a character 
that differs from the one present in mSUGRA-like scenarios. However, it shares similarities 
with gauge mediation, where the gravitino can play a role similar to the neutrino in our 
casen By contrast, the "lightest" neutralino, X^^ ^ typically has non-negligible wino/bino 
components that induce couplings similar to a more standard (massive) neutralino LSP. 
Nevertheless, here this state decays promptly, and is more profitably thought as a neutralino 
LSP in the RPV-MSSM (but with 2-body instead of 3-body decays). 

2.2 The Scalar Sector 

In this section, we discuss the squark, slepton and Higgs sectors, emphasizing the distinctive 
features compared to other supersymmetric scenarios. 

2.2.1 Squarks 

Squarks have interesting non-MSSM properties in the present setup. They are charged 
under the i?-symmetry {R = +1 for the LH squarks and R = —1 for the RH ones), and 
as a result they also carry lepton number. Thus, they are scalar lepto-quarks (strongly 
interacting particles carrying both baryon and lepton number). This character is given by 
the superpotential RPV operator \[-j^LiQjD1, which induces decays such as — )■ hRe\. In 
addition, and unlike in more familiar RPV scenarios, some of these couplings are not free 
but directly related to Yukawa couplings: X'^ = i/d, = Vs and = i/b- The full set of 
constraints on the A' couplings subject to these relations was analyzed in Ref. [l]. The A333 
coupling is the least unconstrained, being subject to 

A'333 < (2.1xlO-2)M ^ 1.4COS/3, (4) 

''There is also a light gravitino in our scenario, but its couplings are suppressed, and plays no role in the 
LHC phenomenology. 
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where yb = "n^^/ve, tan/3 = f^/ve and we took m;,(yU ^ 500 GeV) ~ 2.56 GeV 26 . In 
this work, we will assume that the only non-vanishing A' couplings are those related to the 
Yukawa couplings, together with A333. We will often focus on the case that the upper limit 
in Eq. ^ is saturated, but should keep in mind that A333 could turn out to be smaller, and 
will comment on the relevant dependence when appropriate. 

It is also important to keep in mind that the /^-symmetry forbids any LR mixing. As 
a result, the squark eigenstates coincide with the gauge eigenstates, at least if we neglect 
intergenerational mixing Q We will assume in this work that the first two generation squarks 
are relatively degenerate. As we will see, the current bound on their masses is about 500 — 
700 GeV. We will also see that the third generation squarks can be lighter, possibly consistent 
with estimates based on naturalness from the Higgs sector. 



2.2.2 Sleptons 

The sleptons are expected to be among the lightest sparticles in the new physics spectrum. 
This is due to the intimate connection of the slepton sector with EWSB, together with the 
fact that a good degree of degeneracy between the three generation sleptons is expected. 
The possible exception is the LH third generation slepton doublet, if the RPV coupling A333 
turns out to be sizable. As a result, due to RG running, the LH stau can be several tens of 
GeV lighter than the selectron and smuon, while the latter should have masses within a few 
GeV of each other. Note that the sleptons are -R-neutral, hence do not carry lepton number. 
This is an important distinction compared to the standard extension of lepton number to 
the new physics sector. 

Since the electron sneutrino plays the role of the down-type Higgs, naturalness requires its 
soft mass to be very close to the electroweak scale. To be definite, we take m| ~ m|, ~ (200- 
300 GeV)^. Depending on how this compares to the /x-term, the sleptons can be heavier or 
lighter than the lightest neutralino, . When Xi~ is lighter than the sleptons we will say 
that we have a "neutralino LSP scenario". The other case we will consider is one where the 
LH third generation slepton doublet is lighter than , while the other sleptons are heavier. 
Given the possible mass gap of several tens of GeV between the {v-r-, tl) pair and the other 
sleptons, this is a rather plausible situation. We will call it the "stau LSP scenario" , although 
the r-sneutrino is expected to be up to ten GeV lighter than the stauj^ The possibility that 
several or all the sleptons could be lighter than may also deserve further study, but we 
will not consider such a case in this work. 

We also note that some of the couplings in the RPV operator Xij^LiLjE^ are related to 
lepton Yukawa couplings: A122 = and A133 = y^. The bounds on the remaining Ajj^'s 
under these restriction have been analyzed in jlj, and have been found to be stringent. 

^This assumptions is not necessary, given the mild flavor properties of C/(l)fl-symmetric models [6j|9j[27j. 
This opens up the exciting prospect of observing a non-trivial flavor structure at the LHC, that we leave for 
future work. 

^Again, we remind the reader that we are using standard terminology in a non-standard setting. In 
particular, a rigorous separation of the SM and supersymmetric sectors is not possible, due to the mixings 
in the neutralino and chargino sectors. Also, the supersymmetric particles end up decaying into SM ones, 
similar to RPV-MSSM scenarios. Furthermore, the light gravitino could also be called the LSP, as in gauge- 
mediation. However, unlike in gauge mediation, here the gravitino is very rarely produced in superparticle 
decays, hence not phenomenologically relevant at the LHC. Thus, we will refer to either the [vt^tl) pair or 
as the "LSP" , depending on which one is lighter. Our usage emphasizes the allowed decay modes. 
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We note that, in principle, it could be possible to produce sleptons singly at the LHC 
through the X'-^/^LiQjDl operator, with subsequent decays into leptons via the Xij^LiLjE^ 
induced interactions. We have studied this possibility in Ref [l], and found that there may 
be interesting signals in the ii^ ^~ and e^/i^ channels. However, in this work we do not 
consider such processes any further, and set all A couplings to zero, with the exception of 
the Yukawa ones. The tau Yukawa, in particular, can play an important role. 



2.2.3 The Higgs Sector 

The "Higgs" sector is rather rich in our scenario. The EW symmetry is broken by the vev's 
of the neutral component of the up-type Higgs doublet, and the electron sneutrino, 
Vei which plays a role akin to the neutral down-type Higgs in the MSSM. We have also a 
scalar SM singlet and a scalar SU{2)l triplet, the superpartners of the singlino and tripletino 
discussed in the previous section. These scalars also get non-vanishing expectation values. 
However, it is well known that constraints on the Peskin-Takeuchi T-parameter require the 
triplet vev to be small, f ^ < 2 GeV. We will also assume that the singlet vev is in the few 
GeV range. This means that these two scalars are relatively heavy, and not directly relevant 
to the phenomenology discussed in this paper. Note that all of these states are i?-neutral. 

There is another doublet, i?^, the only state with non-trivial -R-charge (= +2). It does not 
acquire a vev, so that the i?-symmetry is not spontaneously broken, and therefore this state 
does not mix with the previous scalars. Its (complex) neutral and charged components are 
relatively degenerate, with a mass splitting of order 10 GeV, arising from EWSB as well as 
the singlet vev. For simplicity, we will assume its mass to be sufficiently heavy (few hundred 
GeV) that it does not play a role in our discussion. Nevertheless, it would be interesting to 
observe such a state, due to its special -R-charge. 

The upshot is that the light states in the above sector are rather similar to those in the 
MSSM: a light CP-even Higgs, a heavier CP-even Higgs, a CP-odd Higgs and a charged 
Higgs pair. The CP-even and CP-odd states are superpositions of the real and imaginary 
components of and the electron sneutrino (with a small admixture of the singlet and 
neutral triplet states). Given our choice for the slepton soft masses, the heavy CP-even, 
CP-odd and charged Higgses are expected to be relatively degenerate, with a mass in the 
200 — 300 GeV range (the charged Higgs being slightly heavier than the neutral states). The 
charged Higgs is an admixture of and the LH selectron cl (and very suppressed charged 
tripletino components). The RH selectron, as well as the remaining neutral and charged 
sleptons do not mix with the Higgs sector, and can be cleanly mapped into the standard 
slepton / sneutrino terminology. 

The light CP-even Higgs, h, is special, given the observation of a Higgs-like signal by both 



the ATLAS j28] and CMS 29 collaborations at about 125 GeV. This state can also play 
an important role in the decay patterns of the various super-particles. Within our scenario, 
a mass of mh ~ 125 GeV can be obtained from radiative corrections due to the triplet and 
singlet scalars, even if both stops are relatively light (recall the suppression of LR mixing 
due to the i?-symmetry) . This is an interesting distinction from the MSSM. A more detailed 
study of these issue will be dealt with in a separate paper |30j. Here we point out that these 
arguments suggest that Af should be somewhat small, while A^ should be of order one. This 
motivates our specific benchmark choice: Af = and A^ = 1 (although occasionally we will 
allow Af to be non- vanishing) . These couplings affect the neutralino/chargino composition 
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and are therefore relevant for the colhder phenomenology. 



2.3 Summary 

Let us summarize the properties of the superpartner spectrum in our scenario, following from 
the considerations in the previous sections. All the gauginos ( "gluino" , "wino" and "bino" ) 
are relatively heavy, in particular heavier than all the sfermions. The first two generation 
squarks can be below 1 TeV, while the third generation squarks can be in the few hundred 
GeV range. These bounds will be discussed more fully in the remaining of the paper. The 
sleptons, being intimately connected to the Higgs sector, are in the couple hundred GeV 
range. So are the "lightest" neutralino and chargino states, which are Higgsino-like. Mixing 
due to the electron sneutrino vev, induces interesting couplings of the new physics states 
to the electron- neutrino and the electron, while new interactions related to the lepton and 
down-quark Yukawa couplings give rise to non-MSSM signals. The collider phenomenology 
is largely governed by a new (approximately) conserved i?-charge, and will be seen to be 
extremely rich, even though the spectrum of light states does not seem, at first sight, very 
complicated or unconventional. Finally, we mention that there is also an SU{3)c octet scalar 
(partner of the octetinos that are part of the physical gluino states) that will not be studied 



here (for studies of the octet scalar phenomenology, see 31,321). 



3 Sparticle Decay Modes 

In this section we discuss the decay modes of the superparticles relevant for the LHC collider 
phenomenology. We have checked that three-body decays are always negligible and therefore 
we focus on the two-body decays. 



3.1 Neutralino Decays 

From our discussion in the previous section, the lightest (non SM-like) neutralino is a 
Higgsino-like state (that we call X°^), while the truly stable neutralino state is none other 
than the electron-neutrino. It was also emphasized that X^"*" has small, but not always negli- 
gible, gaugino components. The other two (Dirac) neutralino states are heavy. We therefore 
focus here on the decay modes of X^~^. 



As explained in Subsection 2.2.2, we consider two scenarios: a "neutralino LSP scenario", 
where Xi~^ is lighter than the LH third generation slepton doublet, and a "stau LSP scenario" 
with the opposite hierarchy. The decay modes of the lightest neutralino depend on this choice 
and we will consider them separately. 

Neutralino LSP Scenario : 

If X^~^ is lighter than the (i>^, f£) pair, the possible decay modes for have partial decay 
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Figure 3: X^^ branching fractions in the "neutrahno LSP scenario' ' for Mf = 1 TeV Mf = 1.5 TeV, and 
II = 200 GeV. In the left panel we take Af = and = 1, and in the right panel we take Af = Aj = 0.4. 
The former case might be favored by the observation of a Higgs-like state at m/i « 125 GeV. We also take 
the Higgs mixing angles as « 0.98, Rij;, ~ 0.2 and Ris,Rit ^ 1- 



widths [in the notation of Eqs. ([T|-([3])]: 
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(6) 
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2567r 



Xi , 
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(7) 



+ V2 (Af f/il + XlU^,) Vlj^R,^ + V2V,^,Ul {XlRu + X^Ru) 

where we denote the X^^ mass by mj^o, and Ru are the mixing angles characterizing the 
composition of the hghtest Higgs, h. In our scenario all the other Higgs bosons are heavier 
than the lightest neutralino. We note that the above expressions contain an explicit factor of 
l/-\/2 for each occurrence of a neutralino mixing angle, compared to the standard ones 33 



35 



This is because the mixing matrix elements, U,^ and V,'^ are defined in a Dirac basis, 
whereas in the usual approach the neutralinos are intrinsically Majorana particles. Recall 
also that, for simplicity, we are assuming here that all quantities are real. The generalization 
of these and subsequent formulas to the complex case should be straightforward. 

The above decay modes can easily be dominated by the neutrino-neutralino mixing an- 
gles, since the contributions due to the higgsino (U^) and tripletino components are highly 
suppressed. This mixing angles, in turn, are controlled by the sneutrino vev. Note that in 
the RPV-MSSM such decay modes are typically characterized by displaced vertices due to 
the extremely stringent bounds on the sneutrino vev arising from neutrino physics 36 . By 
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contrast, in our scenario the sneutrino vev is allowed to be sizable (tens of GeV), and is in 
fact bounded from below from perturbativity/EWPT arguments, so that these decays are 
prompt. 

The left panel of Fig. |3] shows that the decay width into hug is the dominant one in the 
small sneutrino vev limit, while in the large sneutrino vev limit the channels involving a 
gauge boson can be sizable. We also note that it is possible for the W~e\^ decay channel 
to be the dominant one, as shown in the right panel of Fig. |3j In this case we have chosen 

= = 0.4, which leads to a cancellation between the mixing angles such that Zv^. is 
suppressed compared to W~e\. For such small couplings, the radiative contributions to the 
lightest CP-even Higgs are not large enough to account for the observed rrih ~ 125 GeV, 
while stops (due to the absence of LR mixing) are also not very effective for this purpose. 
Therefore, without additional physics such a situation may be disfavored. We mention it, 
since it is tied to a striking signal, which one should nevertheless keep in mind. 

Stau LSP Scenario: 



If instead the (ut-jT^) pair is lighter than Xi'^, the and Ur^r channels open up with 

partial decay widths given by 
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In Eq. (|8j) we have suppressed addi- 
tional terms proportional to the r Yukawa 
coupling, that give negligible contribu- 
tions compared to the ones displayed. Al- 
though we have included the full expres- 
sions in the numerical analysis, we choose 
to not display such terms to make the 
physics more transparent. The only cases 
where contributions proportional to the 
Yukawa couplings are not negligible occur 
when the top Yukawa is involved]^ We 
then see that Eqs. (|8| and ^ are con- 
trolled by the gaugino components, even 
for the suppressed V-^^ and V^^^^ shown in 
Fig. [2j Thus, these decay channels domi- 
nate over the ones driven by the neutrino- 
neutralino mixing, as shown in Fig. |4j 
Here the v-rVr channel is slightly suppressed compared to the one into the charged lep- 
ton and slepton due to a cancellation between the mixing angles in Eq. ([s]). In other regions 
of parameter space such a cancellation may be more or less severe. 
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Figure 4: X^'^ branching fractions in the "stau LSP 
scenario" for = 1 TeV, Aff = 1.5 TeV, a* = 250 GeV, 
Af = and = 1. We also take m^^ w m^^ = 200 GeV. 
The Higgs mixing angles are as in Fig. |3] 



^Even the contribution from the bottom Yukawa coupling (with possible large tan/3 enhancements) is 
negligible, given the typical mixing angles in the scenario. 
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3.2 Chargino Decays 



The lightest of the charginos (other than the electron) is X^". It is Higgsino-like, which 
follows from its R = —Q nature, and the fact that the winos are heavy. Note that, in 
contrast, the electron and the other charged leptons have R = Q. Therefore, the two- 
body decays of can involve a charged lepton only when accompanied with an electrically 
neutral, \R\ = 2 particle, the only example of which is the R^ scalar. However, this state does 
not couple directly to the leptons We take it to be heavier than X^^, which has important 
consequences for the allowed chargino decay modes. For instance, in the region where is 
heavier than X^^ the potentially allowed decay modes of X^~ are into W^Ue and W^X^~ 
where X^~ denotes the antiparticle of X^'^ . However, the second channel is closed in most of 
the parameter space since X^^ and X^~ are relatively degenerate (with a mass splitting of 
order ten GeV). The dominant decay mode in this "neutralino LSP scenario" has a partial 
decay width given by: 

nxr^W^..) ^ 4(V'r„<-v^niC'5)^m.f (1-^)^(2 + J) (10) 

where we denote the mass of X^^ by m<^±. Therefore, for sufficiently heavy sleptons the 
chargino always decays into W'^Uf.. 

If instead fi is lighter than X^~ one can also have X^ — t- Tj^z/,- with 

r(A-r ^ f>0 = £iU^rra,f (1 - £|j , (11) 

Typically, this decay channel dominates, but the W^Ve can still have an order one branching 
fraction. 



3.3 Slepton Decays 

We focus on the decays of the {pr, tl) pair since it may very well be the "LSP", i.e. the last 
step in a cascade decade to SM particles. In this case the charged slepton decay modes are 
— )■ T^Ve and — )• lihR, with partial decay widths given by: 

nfi^v.) ^'^vl' (12) 
mr^f.M - ?^(AW^(l-^)^ (13) 

The decay widths for the SU{2)l related processes, i>r — ?• T/jC^ and Vr b^bji, are obtained 
from Eqs. (12) and (13) with the replacements nif^ — )■ rrij^^ and rrit — )■ m^. In Fig. [sj we 



show the branching fractions as a function of the sneutrino vev, assuming that A333 saturates 
Eq. ([4]), and taking rrir = 1.7 GeV. We see that the tLbji channel can be sizable in the large 
sneutrino vev/small tan /3 limit, in spite of the phase space suppression when m^^ ~ nit + mi, 
(left panel). Away from threshold, it can easily dominate (right panel). 



''Recall that the Rd SU{2) doublet does not play any role in EWSB. 
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Figure 5: f^, (solid lines) and Vt (dashed lines) decay modes for two masses: m^f^ w m^^ — 180 GeV (left 
panel) and ra^^ rat,^ — 250 GeV (right panel). It is assumed that is heavier than the {vt^t'l) pair, 
and that A333 saturates Eq. Q. 



If, on the other hand, and X^^ are hghter than the LH third generation sleptons, 
their dominant decay modes would be fj^ — >■ X^t]] or — )■ X\ for the charged lepton, 
and Vt — )■ X\^Vt or — )■ X^^r^ for the sneutrino. 



3.4 Squark Decays 

As aheady explained, we focus on the case where the gluinos are heavier than the squarks 
and, therefore, the squark decay mode into a gluino plus jet is kinematically closed. The 
lightest neutralinos and charginos are instead expected to be lighter than the squarks since 
naturalness requires the /x-term to be at the electroweak scale, while we will see that the first 
and second generation squarks have to be heavier than about 600 GeV. Thus, the squark 
decays into a quark plus the lightest neutralino or into a quark plus the lightest chargino 
should be kinematically open. However, the decay mode of the left handed up-type squarks, 
which have Q = 2/3 and i? = 1, into the lightest chargino X~l plus a (i?-neutral) jet is 
forbidden by the combined conservation of the electric and R-charges: ui^ -j^ X^~ j. The 
decay mode into the second lightest neutralino, which can be of the (++) type, could be 
allowed by the quantum numbers, but our choice > rriq ensures that it is kinematically 
closed. Note also that since ur has Q = 2/3 and R = —1, one can have ur — t- X^^j. 



3.4.1 First and Second Generation Squarks 

• The left-handed up-type squarks, ul and cl, decay into X^^j and e^j with: 
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(14) 
(15) 



and analogous expressions for cl (in Eq. (14), we do not display subleading terms 
proportional to the Yukawa couplings). The second decay is an example of a lepto- 
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quark decay mode. However, taking into account the smallness of the Yukawa couphngs 
for the first two generations, together with the Xi~^ composition shown in Fig. ^ one 
finds that the dominant decay mode is the one into neutrahno and a jet. Therefore, 
in the region of parameter space we are interest in, ul and cl decay into X^~^j with 
almost 100% probabihty. 



The down-type left-handed squarks, di and Sl, have the following decay channels: 
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with analogous expressions for sl- The relative minus sign in the gaugino contributions 
to the neutralino decay channel is due to the SU{2) charge of the down-type squarks. 



and should be compared to the up-type case, Eq. (14). This leads to a certain degree 



of cancellation between the contributions from the bino and wino components, which 
together with the factor of 1/18 results in a significant suppression of the neutralino 
channel. Since the Yukawa couplings are very small, it follows that the chargino channel 
is the dominant decay mode of the down-type squarks of the first two generations. 



The right-handed up-type squarks, ur and cr, decay according to 
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(19) 
(20) 



with analogous expressions for cr. The chargino decay mode fo ur is suppressed since 
the up-type Yukawa coupling is very small. Therefore, the right-handed up-type squark 
decays into X^^ j with almost 100% probability. However, the charm Yukawa coupling 



is such that the various terms in Eqs. (19) and (20) are comparable when the mixing 



angles are as in Figs. [T] and |2j For this benchmark scenario, both decay channels 
happen to be comparable, as illustrated in the left panel of Fig. [6j Here we used 
'-vl with mc{^l ^ 600 GeV) ^ 550 MeV (261. 
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The right-handed down-type squarks, dR and sr, decay according to 
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Figure 6: Branching fractions for cr (left panel) and sr (right panel) taking Aff = 1 TcV, = 1.5 TeV, 
= 200 GeV, Af = and = 1. 
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Figure 7: Branching fractions for the decay modes computed for Mf 1 TeV, = 1.5 TeV, 

fi ~ 200 GeV, Af = and Aj = 1. We also assume A333 = (2.1 x 10^^)/?/b. In the left panel we take 
rrij:^ — 500 GeV, and show the dependence on the sneutrino vev. In the left panel we show the dependence 
on rrij:^ for Ve = 10 GeV (solid lines) and Vg = 50 GeV (dashed lines). 



with analogous expressions for Sr. Again, for the down squark the Yukawa couphngs 
are neghgible so that it decays dominantly into neutrahno plus jet. For the strange 
squark, however, the various channels can be competitive as illustrated in the right 



panel of Fig. M Here we used i/s = ms/vg with ms{n f» 600 GeV) ^ 49 MeV 26 



3.4.2 Third Generation Squarks 

For the third generation we expect the lepto-quark signals to be visible in all of our parameter 
space, although they may be of different types. The point is that the bottom Yukawa coupling 
can be sizable in the small sneutrino vev/large tan/3 limit (as in the MSSM), thus leading 
to a signal involving first generation leptons through the A'j^gg = yb ^ 1.15 x 10~^sec/3 
coupling. In the large sneutrino vev/small tan /3 limit, on the other hand, the RPV coupling 
A333 < 1.4 cos /3 can be of order of g', and may lead to third generation leptons in the final 
state. 
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The left-handed stop, ti, has the following decay modes: 
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When kinematically allowed, the decay mode into neutralino plus top is the dominant 
one since it is driven by the top Yukawa coupling, as shown in Fig. [7j However, this 
figure also shows that the two lepto-quark decay modes can have sizable branching 
fractions]^ In particular, at small sneutrino vev the electron-bottom channel is the 
dominant lepto-quark decay mode (since it is proportional to the bottom Yukawa), 
while in the large vev limit the third generation lepto-quark channel dominates [we 
have assumed that A333 saturates the upper bound in Eq. Q]. The existence of lepto- 
quark channels with a sizable (but somewhat smaller than one) branching fraction 
is a distinctive feature of our model, as will be discussed in more detail in the fol- 
lowing section. We also note that in the case that A333 is negligible and does not 
saturate the bound in Eq. the ii ^ Tl^r channel is no longer present, so that the 
BR(ti ^ ^l^r) c'-iid BR(ti — )■ X'l'^t) increase in the large sneutrino vev limit (but are 
qualitatively the same as the left panel of Fig. [T]). 

The left-handed sbottom, b^, has several decay modes as follows: 
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When kinematically open, the dominant decay mode is into a chargino plus top since 



*Here we used yt = mt/^/v'^ — and j/b = rni,/vg with rritifJ. ~ 500 GeV) « 157 GeV and mf,(/i 
500 GeV) « 2.56 GeV [26]. 



17 




Figure 8: Branching fractions for 6l computed for A/f = 1 TeV, = 1.5 TeV, ^l = 200 GeV, Af = 
and = 1. We also take A333 = (2.1 x 10^-^)/?/f,, and add together the two neutrino channels {De and i>r)- 
In the left panel we take m-j^^ = 500 GeV, and show the dependence on the sneutrino vev. In the left panel 
we show the dependence on mj,^ for Ve = 10 GeV (solid lines) and Vg = 50 GeV (dashed lines). 



it is controlled by the top Yukawa coupling. The decays into neutrino plus bottom 
have always a sizable branching fraction, as can be seen in Fig. |8| However, one should 
note that when A333 is negligible, so that the 6l — )■ i^rbn channel is unavailable, the 



decay involving a neutrino {u^ 
order 0.3% at Vg = 50 GeV). 
change qualitatively. 



only) decreases as the sneutrino vev increases (being of 
The other two channels adjust accordingly, but do not 



For the right-handed stop, tji, the decay widths are: 
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For the benchmark choice of Mf = 1.5 TeV, Mf = 1 TeV, 12 = 200 GeV, Af = and 
= 1, we have r(t|j ^ ^1°+^^) = 26% (15%) and r(t^ ^ X^-Br) = 74% (85%) for 
m^^ = 500 (400) GeV, independently of the sneutrino vev. For m^^ < m^o + nit, the 

RH stop decays into X^~bfi essentially 100% of the time. See left panel of Fig. lo] 
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Figure 9: Branching fractions for as a function of m^-^ (left panel), and for as a function of Vg (right 
panel) computed for Aff = 1 TeV, = 1.5 TeV, = 200 GeV, Af = and = 1. For Ir, we take 
A333 = (2.1 X 10^^)/?/h, assume m^^^ ^ m^o,mt, and add together the two neutrino channels (j>e and ;>,-)■ 



The right-handed sbottom, b^, has a variety of decay modes: 
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nh-^r-t,) = (^333^ (^1 - ^ j , (37) 

nh^UrbL) = ^^(^333)'- (38) 

The lepto-quark signals are the dominant ones. Adding the two neutrino channels, the 
decay mode into ub has a branching fraction of about 50% as shown in the right panel 
of Fig. |9j The charged lepton signals can involve a LH electron or a r plus a top quark. 
Note also that the decay mode into X^^b is very suppressed. We finally comment on 
the modifications when A333 is negligible. Once the b^ — )■ r^ti and bn — J- i^rbi channels 
become unavailable, one has that BR{bji — )■ ej^ti) ~ 0.6 and BR(6ij — )• Uebi) ~ 0.4, 
independent of the sneutrino vev. The b*j, — X^^bji channel remains negligible. 



4 and 2^^ Generation Squark Phenomenology 

In the present section we discuss the LHC phenomenology of the first and second generation 
squarks, which are expected to be the most copiously produced new physics particles. Al- 
though these squarks are not required by naturalness to be light, flavor considerations may 
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Figure 10: Left Panel: qq* production cross-section (all flavor combinations), for the 7 TeV LHC run, 
computed for 2 TeV Dirac (red) and Majorana (blue) gluinos. In the right panel we plot the ratio between 
the two cross sections, showing that the suppression in the Dirac case can be significant. 



suggest that they should not be much heavier than the third generation squarks. Therefore, 
it is interesting to understand how hght these particles could be in our scenario. As we will 
see, current bounds allow them to be as light as 500 — 700 GeV, while in the MSSM the 
LHC bounds have already exceeded the 1 TeV threshold. The bounds can arise from generic 
jets + searches, as well as from searches involving leptons in the final state. 



4.1 Squark Production 



We compute the cross section to produce a given final state X in our model as follows: 

a{pp ^ X) = a{pp ^i)x BR(i ^ X) , (39) 

i 

where i = qiq2, gq, gg, and the squark pair production can in principle come in several flavor 
and chirality combinations. We generate the production cross section for each independent 
i-th state with MadGraphS 37 . Here we note that, due to the assumption of gluinos 
in the multi-TeV range, and the fact that we will be interested in squarks below 1 TeV, 
our cross section is dominated by the production of squark pairs. We have also computed 

as a function of the squark mass for fixed 
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the corresponding ii'-factor with Prospino2 
(Majorana) gluino masses of 2 — 5 TeV. We find that for squark masses below about 1 TeV, 
the K-factor is approximately constant with K ^ 1.6. Since, to our knowledge, a NLO 
computation in the Dirac case is not available, we will use the previous K-factor to obtain a 
reasonable estimate of the Dirac NLO squark pair-production cross-section. 

One should note that the Dirac nature of the gluinos results in a significant suppression 
of certain t-channel mediated gluino diagrams compared to the Majorana (MSSM) case, 
as already emphasized in 14 , T5] (see also Fig. 10). Nevertheless, at = 2 TeV such 



15] (see also Fig. 

contributions are not always negligible, and should be included. For instance, we find that 
for degenerate squark with nig = 800 GeV, the production of ulUr, Uidji and URdi is 
comparable to the "diagonal" production oidjffl and qRq*R for all the squark flavors q = 
u, d, s, c taken together. As indicated in Eq. (39) we include separately the BR for each i-th 



state to produce the final state X, since these can depend on the squark flavor, chirality or 
generation. 
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4.2 "Simplified Model" Philosophy 

We have seen that ul, uji, da and cl decay dominantly through the neutrahno channel, the 
LH down- type squarks, di and sl^ decay dominantly through the chargino channel, and cr 
and sb. can have more complicated decay patterns (see Fig. [6]). The striking lepto-quark decay 
mode, sr — )■ e^j, will be treated separately. In this section we focus on the decays involving 
neutralinos and charginos. Since the signals depend on how the neutralino/chargino decays, 
it is useful to present first an analysis based on the simplified model (SMS) philosophy. To 
be more precise, we set bounds assuming that the neutralinos/charginos produced in squark 
decays have a single decay mode with BR = 1. We also separate the "neutralino LSP 
scenario", in which X^'^ /X^~ decay into SM particles, from the "stau LSP scenario", where 
they decay into f^r^, Uri^r or f^i'r- We will give further details on these subsequent decays 
below, where we treat the two cases separately. 

Here we emphasize that we regard 
the jets plus X^'^ /X^" stage as part of 
the production. The point is that an im- 
portant characteristic of our scenario is 
that different types of squarks produce 
overwhelmingly only one of these two 
states. For instance, if we are interested 
in two charginos in the squark cascade 
decays, this means that they must have 
been produced through LH down-type 
squarks (with a smaller contribution 
from CrC*j^ production), and the produc- 
tion of any of the other squarks would 
not be relevant to this topology. Con- 
versely, if we are interested in a topol- 
ogy with two neutralinos, the LH down- 
type squarks do not contribute. We de- 
note by (Ti the corresponding cross sec- 
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tions, computed via Eq. (39) with X = 
"X+-Xr+jj" or X = "XrX?-jj", 
taking the BR's as exactly zero or one, 
according to the type of squark pair 
At the same time, since in other realiza- 
tions of the i?-symmetry these produc- 



Figure 11: Cross-sections for the separate production 
of X°+X°-, + and X°+X+- trough squark pair- 

production in the SMS approach (see main text). The solid 
and dotted hnes correspond to cti, according to the case. 
The dashed hne marked as (72 corresponds to the full pair- 
production of squarks, irrespective of how they decay. The 
cross-sections are computed for Mg = 2 TeV for a 7 TeV 
LHC run, with a K-factor, K = 1.6. 



tion patterns may not be as clear-cut, 

we will also quote bounds based on a second production cross section, denoted by a2, where 
it is assumed that all the squarks decay either into the lightest neutralino or chargino chan- 
nels with unit probability. This second treatment is closer to the pure SMS philosophy, 
but could be misleading in the case that lepton number is an i?-symmetry. We show the 



corresponding cross-sections in Fig. 11 



X 



^Thc only exception is the RH charm squark, cji, for which we take BR(c|j -^^i^j) = BR(cj^ — > 
^ j) = 0.5, although the characteristics of the signal are not very sensitive to this choice. We also neglect 



the decays of Sf; into neutralino/ncutrino plus jet. 
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It should also be noted that the great majority of simplified models studied by ATLAS 
and CMS consider either nig = Mg, or niq ^ Mg. Therefore, at the moment there are only a 
handful of dedicated studies of our topologies, although we will adapt studies performed for 
other scenarios to our case. In the most constraining cases, we will estimate the acceptance by 
simulating the signal in our scenaricp°1 and applying the experimental cuts, but for the most 
part a proper mapping of the kinematic variables should suffice (provided the topologies 
are sufficiently similar). A typical SMS analysis yields colored-coded plots for the upper 
bound on (T X BR (or Axe) for the given process, in the plane of the produced (strongly- 
interacting) particle mass (call it rriq), and the LSP mass (call it mLSp)- Iii most cases, the 
LSP is assumed to carry ^t- Often, there is one intermediate particle in the decay chain. Its 
mass is parametrized in terms of a variable x defined by ^intermediate = xrrig + (1 — x)mLSP- 
In our "neutralino LSP scenario", the intermediate particle is either the lightest neutralino 
Xi~^ or the lightest chargino , whose masses are set by the //-term. Since the particle 
carrying the is the neutrino, i.e. mLSP = 0, we have x ~ fi/rrig. 

We will set our bounds as follows: we compute our theoretical cross section as described 
above (i.e. based on the ai or o"2 production cross-sections) as a function of the squark 
mass, and considering the appropriate decay channel for the X^'^ / X^^ (with BR = 1 in 
the SMS approach). Provided the topology is sufficiently similar, we identify the x-axis on 
the color-coded plots in the experimental analyses (usually nig) with nig, take m^sp = 
(for the neutrino), and identify "x" as n/nig (from our discussion above). Then, we increase 
the squark mass until the theoretical cross-section matches the experimental upper bound, 
defining a lower bound on niq. In a few cases that have the potential of setting strong 
bounds, but where the experimentally analyzed topologies do not exactly match the one in 
our model, we obtain the signal e x A from our own simulation and use the 95% CL. upper 
bound on 0" X e X A to obtain an upper bound on a that can be compared to our model 
cross-section. If there are several signal regions, we use the most constraining one. 



4.3 Neutralino LSP Scenario 

In the neutralino LSP scenario, and depending on the region of parameter space (e.g. the 
sneutrino vev or the values of the Xs and At couplings), the lightest neutralino, Xi~^, can 
dominantly decay into Zug, hv^ or W~e\. The "lightest" chargino X^ always decays 
into W^Ue- Following the philosophy explained in the previous subsection, we set separate 
bounds on four simplified model scenarios: 

(1) q ^ (Zz/e) 3 , (3) q ^ X?+j ^ (W^'e^) j , 

(2) q ^ ^ {hu,) J , (4) q ^ X+'J ^ (W+u,) j , 



as well as on two benchmark scenarios (to be discussed in Subsection 4.3.1) that illustrate 
the bounds on the full model. 

There are several existing searches that can potentially constrain the model: 



^°We have implemented the fuh model in FeynRules 38 , which was then used to generate MadGraph 5 
code |37|. The parton level processes are then passed through Pythia for hadronization and showering, and 
through Delphes 



39 



for fast detector simulation. 
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Table 1: Bounds on 1*** and 2"*^ generation squark masses from squark pair production in the "neutralino 
LSP scenario" for the SimpUfied Models (l)-(4), and two benchmark scenarios. See text for further details. 



• jets + , • OS dileptons + + jets , 

• 1 lepton + jets + , • multilepton + jets + (with or without Z veto). 

• Z{U) + jets + , 

We postpone the detailed description of how we obtain the corresponding bounds to 
the appendix, and comment here only on the results and salient features. We find that 
typically the most constraining searches are the generic jets + searches, in particular the 



most recent ATLAS search with 5.8 fb~ 19 . In addition, some of the simplified topologies 
can also be constrained by searches involving leptons + jets + ^t- For example, those 
involving a leptonically decaying Z are important for the Xi~^ — )■ Zug case, while a number 
of multi-lepton searches can be relevant for the topologies that involve a W. We summarize 
our findings in Table [T| where we exhibit the searches that have some sensitivity for the 
given SMS topology. We show the lower bounds on the squark masses based on both the 



0"! and (72 production cross-sections, as described in Subsection |4.2[ We see that these are 
below 650 GeV (based on ai; the bound from (T2 is provided only for possible application 
to other models). We also show the bounds for two benchmark scenarios (which depend 
on the sneutrino vev), as will be discussed in the next subsection. These are shown under 
the (Ji column, but should be understood to include the details of the branching fractions 
and various contributing processes. We have obtained the above results by implementing 
the experimental analysis and computing the relevant e x A from our own simulation of 
the signal, and using the model-independent 95% CL upper bounds on cr x e x A provided 
by the experimental analysis. Whenever possible, we have also checked against similar 
simplified model interpretations provided by the experimental collaborations. Such details 
are described in the appendix, where we also discuss other searches that turn out to not 
be sensitive enough, and the reasons for such an outcome. In many cases, it should be 
possible to optimize the set of cuts (within the existing strategies) to attain some sensitivity. 
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Figure 12: Cross-sections for a variety of signatures in the "neutralino LSP scenario". All are computed 
for = 1 TeV, = 1.5 TeV, ^l = 200 GeV, and assuming that A333 saturates Eq. Q. In the left panel 
we take niq = 700 GeV and Af = 0, = 1 (benchmark 1), while in the right panel we use = 550 GeV 
and Af = A„ = 0.4 (benchmark 2). 



This might be interesting, for example, in the cases involving a Higgs, given that one might 
attempt to reconstruct the Higgs mass. 

We turn next to the analysis of the full model in the context of two benchmark scenarios. 

4.3.1 Realistic Benchmark Points 

Besides the "simplified model" type of bounds discussed above, it is also interesting to 
present the bounds within benchmark scenarios that reflect the expected branching fractions 



for the neutralinos/charginos discussed in Subsections 3.1 and 3.2 One difference with the 
analysis of the previous subsections is that we can have all the combinations of X^~^ X^^ j j , 
X^^ X^^jj and X^^X^^jj in squark decays, with the corresponding BR's. In Fig. |ll| we 
have shown the individual cross-sections in the SMS approach. These give a sense of the 
relative contributions of the various channels. In particular, we see that the channel 
dominates. 

Benchmark 1 : (Mf =^1 TeV Mf = 1.5 TeV, = 200 GeV, Af = 0, = 1) corresponds 
to the case that the X^'^ — )■ hve decay channel is important (in fact, dominant at small 
sneutrino vev), while the gauge decay channels of the X^'^ can be sizable (see left panel of 
Fig. [3]). The LHC searches relevant to this scenario are: 

• jets + 5 • dileptons (from Z decay) + jets + , 

• 1 lepton + jets + , • multilepton + jets + (without Z cut). 

• OS dileptons + + jets , 
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We apply the model-independent bounds discussed in the previous sections, and find 
that the jets + search is the most constraining one. Using o^j^Xt ~ 20 — 40 fb, we find 
mg- > 620 - 690 GeV [m-^ > 590 - 650 GeV) for Vg = 10 GeV (vg = 50 GeV). We show in the 
left panel of Fig. IT2| the cross-sections for several processes, for rriq = 700 GeV. These are 



computed from Eq. (39) using the actual BR's for the chosen benchmark. Although there is 
some dependence on the sneutrino vev, the global picture is robust against Vg. 

Benchmark 2 : (Mf = 1 TeV, Mf = 1.5 TeV, /x = 200 GeV, = = 0.4) corresponds 
to the case that the Xj'^ — W"e^ decay channel dominates (see right panel of Fig. p). 



In the right panel of Fig. 12, we show the cross-sections for the main processes. We see 
that, for this benchmark, the "leptonic channels" have the largest cross sections (especially 
the multilepton -|- jets + one). However, taking into account efficiencies of at most a 
few percent for the leptonic searches (as we have illustrated in the previous section), we 
conclude that the strongest bound on the squark masses arises instead from the jets + 
searches (as for benchmark 1). Using ct^+Xt ~ 20 — 40 fb, we find nig > 520 — 580 GeV 
{rrig > 500 - 560 GeV) for Ve = 10 GeV (ve = 50 GeV). Note that there is a sizable "no 
missing energy" cross section. However, this could be significantly lower once appropriate 
trigger requirements are imposed. 

4.4 Stau LSP Scenario 

In this scenario the dominant decay modes of X^'^ are into r^r^ or Uri^r (about 50-50), 
while the chargino X^' decays into t/z^t- The decay modes of f£ depend on the sneutrino 
vev: for large Vg it decays dominantly into tifei? (assuming A333 is sizable), while for smaller 
values of Vg it decays dominantly into r^Uf. trough the r Yukawa coupling. Similarly, z/^ 
decays into 6l6_r for large sneutrino vev, and into T^e^ for small sneutrino vev. In the "stau 
LSP scenario" we prefer to discuss the two limiting cases of small and large sneutrino vev. 



rather than present SMS bounds (recall from Fig. 11 that the squarks produce dominantly 
X°+X°~ pairs). This scenario is, therefore, characterized by third generation signals. 

4.4.1 f£ — )■ T^Ue and i)r Tr&\ decay modes 

These decays are characteristic of the small sneutrino vev limit. In this case all the final 
states would contain at least two taus: i) for the X^'^X^~ topology the final state contains 
2 jets, missing energy and 2r + 2e, 3r + le or 4r's; ii) for the X^'^X^~ topology the final 
state contains 2 jets, missing energy and 2r + le or 3r's; Hi) for the X^^Xi~ topology the 
final state contains 2 jets, missing energy and 2r's. It is important that cases i) and ii) can 
be accompanied by one or two electrons, given that many searches for topologies involving 
r's impose a lepton (e or /i) veto. 



Thus, for instance, a recent ATLAS study 44 with 4.7 fb~^ searches for jets + 
accompanied by exactly one (hadronically decaying) r + one lepton (e or /x), or by two 
r's with a lepton veto. Only the former would apply to our scenario, setting a bound of 



cr X e X A = 0.68 fb"^ A previous ATLAS search 45 with 2.05 fb ^ searches for at least 
2r's (with a lepton veto), setting a bound oi a x e x A = 2.9 fb~^. However, the efficiency 
of such searches is lower than the one for jets plus missing energy (also with lepton veto). 



Understood as hadronic r's. 
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Since in our scenario the cross sections for these two signatures is the same, the latter will 
set the relevant current bound. 



There is also a CMS study 46 sensitive to 4r signals in the context of GMSB scenarios, 
which has a similar topology to our case (SMS: gg production with g — )■ qqXi and Xi 
T^T^Gfj,). From their Fig. 9b, we can see that the 95% CL upper limit on the model cross 
section varies between 0.3 — 0.03 pb for 400 GeV < rrig < 700 GeV. Including the branching 



fractions, and reinterpreting the bound in the squark mass plane, ^ we find a bound of 
rriq > 600 GeV at Vg = 10, where the cross section is about 45 fb. When the sneutrino vev 
increases the bound gets relaxed so that for Vg > 20 GeV there is no bound from this study. 

The generic searches discussed in previous sections (not necessarily designed for sensitiv- 
ity to the third generation) may also be relevant: 

• jets + Xt , • jets + + OS dileptons , 

• jets + + 1 lepton , • jets + + multi leptons , 

• jets + + SS dileptons , 

where the leptons may arise from the Uj. decay as in cases i) and ii) above, or from leptonically 



decaying r'sj^J^ It turns out that, as in the "neutralino LSP scenario", the strongest bound 
arises from the jets + search. We find from simulation of the signal efficiency times 



acceptance for the ATLAS analysis 19 in our model that the most stringent bound arises 



from signal region C (tight), and gives an upper bound on the model cross section of about 
120 fb. Thus, we find that rUg > 500 GeV for Vg = 10 GeV. 



4.4.2 — )• tibji and i>T hLbR decay modes 

When the third generation sleptons decay through these channels, as is typical of the large 
sneutrino vev limit, the signals contain a hh and/or a ti pair, as well as r's. Note that when 
the r's and tops decay hadronically one has a signal without missing energy. However, the 
branching fraction for such a process is of order BR(g — > X°^j)^ x BR(Xj'^ — )► f^r^Y x 
BR(f^ ^ tLbR)^ X BR(t ^ bW+y x BR{W jjf x BR(r ^ jjf ~ few per cent (in the 
large sneutrino vev limit, where all of these branching fractions are sizable). Indeed, we find 
that the "no Xr" cross section for 700 GeV squarks in the "stau LSP scenario" is of order 
1 fb, which is relatively small. Rather, the bulk of the cross section shows in the jets + 
and 1 lepton + jets + channels (with a smaller 2 lepton + jets + contribution). 
Simulation of the ATLAS ]+^t search [l9] in this region of our model indicates that again 
the most stringent bound arises from signal region C (tight) of this study, and gives an 
upper bound on the model cross section of about 70 fb. This translates into a bound of 
niq > 550 GeV for Ve = 50 GeV. 

^^As usual, the topology of this study contains two additional hard jets at the parton level compared to 
our case. 

^■^Note that when there are two taus and no additional electrons, the SS dilepton searches do not apply. 
This is a consequence of the conserved i?-symmetry. 
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5 Third Generation Squark Phenomenology 



We turn now to the LHC phenomenology of the third generation squarks. We start by 
studying the current constraints and then we will explain how the third generation provides 
a possible smoking gun for our model. We separate our discussion into the signals arising 
from the lepto-quark decay channels, and those that arise from the decays of the third 
generation squarks into states containing X^^ or (or their antiparticles). 



5.1 Lepto-quark Signatures 

Due to the identification of lepton number as an i?-symmetry, there exist lepto-quark (LQ) 
decays proceeding through the LQD^ couplings. These can be especially significant for 
the third generation squarks. As discussed in Subsection |3.4.2 in our scenario we expect: 



tL e'lbR, ti Tlbn, bi (z^e + z^t)&r, bn (z/g + e^ti and bn ti- It 



may be feasible to use the channels involving a top quark in the final state 47 , but such 



searches have not yet been performed by the LHC collaborations. Thus, we focus on the 



existing eejj 48,49 , uubb 50 and rrbb [51] searches, where in our case the jets are really 
6-jets|^ The first and third searches have been performed with close to 5 fb~^ by CMS, while 



the second has been done with 1.8 fb~ . In the left panel of Fig. [13| we show the bounds 
from these searches on the LQ mass as a function of the branching fraction of the LQ into 
the given channel. The bounds are based on the NLO strong pair-production cross-section. 
We see that the most sensitive is the one involving electrons, while the one involving missing 
energy is the least sensitive. This is in part due to the lower luminosity, but also because in 
the latter case the search strategy is different since one cannot reconstruct the LQ mass. 



In the right panel of Fig. 13 we show the corresponding branching fractions in our scenario 
as a function of the sneutrino vev, assuming ttilq = 400 GeV (which, as we will see, turns 
out to be the mass scale of interest). We have fixed = 1.5 TeV and Mf^ = 1 TeV, 
and scanned over fi e [—300, 300] GeV and , G [0, 1], which is reflected in the width of 
the bands of different colors. We assume that A333 saturates Eq. (|4]). The BR's are rather 
insensitive to Af and A^, but depend strongly on fi, especially when > 200 GeV. The 
reason is that for larger n the neutralinos and charginos become too heavy, the corresponding 
channels close, and the LQ channels can dominate. This affects the decays of tL and bi, but 
not those of as can be understood by inspecting Figs. ^ [s] and the right panel of Fig. of^ 



The darker areas correspond to the region \n\ G [0,200], while the lighter ones correspond 
to G [200,300]. We can draw a couple of general conclusions: 

1. The uubb branching fractions are below the sensitivity of the present search, except 
when the neutralino/chargino channels are suppressed or closed for kinematic reasons. 
Even in such cases, the lower bound on m^^ is at most 350 GeV. Note that bR is 
unconstrained. 



^^It would be interesting to perform the eejj search imposing a 6-tag requirement that would be sensitive 
to our specific signature. 

^^Note, in particular, that the neutralino decay channel of bji is always suppressed, so that its branching 
fractions are insensitive to fi, unlike the cases of II and b^. This is why the "6jj — (i^g + J^r)&L band" in 



Fig. 13 appears essentially as a line, the corresponding BR being almost independent of /i. 
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Figure 13: Left panel: current bounds on lepto-quark masses ft-om three channels: eejj (blue), rrbb (red) 
and wbh (green) , as a function of the lepto-quark branching fraction into the corresponding channel (based on 
the CMS analyses 49 -51 ). Right panel: Branching fractions into lepto-quark channels for ttilq = 400 GeV, 
as a function of the sneutrino vev, for Mf — 1 TeV, = 1.5 TeV, and scanning over Af, e [0,1] 



and € [0,200] GeV (darker areas) or 
involving a top quark. 



/i| G [200,300] GeV (lighter areas). We do not show the channels 



2. The rrbb search, which is sensitive to BR's above 0.3, could set some bounds at large 
Vg in some regions of parameter space. Such bounds could be as large as 520 GeV, but 
there is a large region of parameter space that remains unconstrained. 

3. The eejj search , which is sensitive to BR's above 0.15, could set some bounds at small 
Vg in some regions of parameter space. Such bounds could be as large as 815 GeV if 
Vg ~ 10 GeV and the neutralino/chargino channels are kinematically closed. However, 
in the more typical region with /i < 200 GeV the bounds reach only up to 550 GeV 
in the small Vg region. Nevertheless, there is a large region of parameter space that 
remains completely unconstrained. 

The latter two cases are particularly interesting since the signals arise from the (LH) 
stop, which can be expected to be light based on naturalness considerations. In addition to 
the lessons from the above plots, we also give the bounds for our benchmark scenario with 
Mf = 1 TeV, Mf = 1.5 TeV, /i = 200 GeV, Af = and = 1, assuming again that Eq. (g) 
is saturated. We find that the uubb search requires m^^ > 350 GeV, and gives no bound on 
TTL-j^^- The rrbb search gives a bound on m^^ that varies from 380 to 400 GeV as Vg varies from 
20 — 50 GeV. The eejj search gives a bound on m^^ that varies from 470 down to 300 GeV 
as Vg varies from 10 — 30 GeV. The other regions in Vg remain unconstrained at present. In 
our benchmark, when m^^ ~ 350 GeV, we expect m^^ ~ 380 — 390 GeV, depending on the 
scalar singlet and (small) triplet Higgs vevs (and with only a mild dependence on Vg). We 
conclude that in the benchmark scenario a 400 GeV LH stop is consistent with LQ searches. 
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while offering the prospect of a LQ signal in the near future, possibly in more than one 
channel. 

Comment on LQ signals from 2"^*^ generation squarks 

We have seen that the RH strange squark has a sizable branching fraction into the LQ 



channel, s/j — )■ e^j, of order 0.4 — 0.65. From the left panel of Fig. 13, we see that the eejj 
CMS lepto-quark search gives a bound of rris ~ 530 — 630 GeV, which is quite comparable to 
(but somewhat weaker than) the bounds obtained in Section |4| Thus, a LQ signal associated 
to the RH strange squark is also an exciting prospect within our scenario. 



5.2 Other Searches 

There are a number of searches specifically optimized for third generation squarks. In addi- 
tion, there are somewhat more generic studies with b-tagged jets (with or without leptons) 
that can have sensitivity to our signals. We discuss these in turn. 

Direct stop searches: In the case of the top squark, different strategies are used to sup- 
press the ti background depending on the stop mass. However, the searches are tailored to 
specific assumptions that are not necessarily satisfied in our framework: 



Perhaps t 
search 
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le most directly applicable search to our scenario is an ATLAS GMSB 
(iiil pair production with ii — )■ txi or ii — )■ bxt and finally Xi ~^ ZG 
or Xi ~^ W^G), so that the topologies are identical to those for LH and RH stop 
pair-production in our model, respectively, with the replacement of the light gravitino 
by Ue (although the various branching fractions are different; see Figs. [7]and|9]for our 



benchmark scenario). Ref. 52 focuses on the decays involving a Z, setting bounds of 
0" X e X A = 18.2 (9.7) fb for their signal region SRI (SR2). Simulation of our signal for 
our benchmark parameters and taking 400 GeV LH stops gives e ^ A 1.9% (1.7%) 
for SRI (SR2), which include all the relevant branching fractions. The corresponding 
bound on the model cross section would then be ci^i^ ~ 1 (0.6) pb. However, a cross- 
section of 0.6 pb is only attained for stops as light as 300 GeV, and in this case the 
efficiency of the search is significantly smaller, as the phase space for the ti X°^t 
decay closes (recall that due to the LEP bound on the chargino, and the Higgsino-like 
nature of our neutralino, the mass of X^^ must be larger than about 100 GeV). We 
conclude that this search is not sufficiently sensitive to constrain the LH stop mass. 
Also, the requirement that the topology contain a Z gauge boson makes this search 
very inefficient for the RH stop topology: t/j — ?■ X^^bji, X^^ — )■ W^Vf., so that no 
useful bound can be derived on m^. 

There is a search targeted for stops lighter than the top {t — )■ hxt , followed by Xi ~^ 
W'^xV)- This is exactly the topology for production in our scenario (with mLSP = 
for the neutrino), but does not apply to tz, since its decays are dominated by lepto- 



quark modes in this mass region. Fig. 4c in 53 shows that for a chargino mass of 
106 GeV, stop masses between 120 and 164 GeV are excluded. As the chargino mass 
is increased, the search sensitivity decreases, but obtaining the stop mass limits would 
require detailed simulation in order to compare to their upper bound cr x e x A = 
5.2 - 11 fb. 
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There are also searches for stop pair production with t — )■ txi- A search where both 



tops decay leptonically 54 would yield the same final state as for t^t^ production in 
our case {bbW^W~ with the W^s decaying leptonically). However, the kinematics 
is somewhat different than the one assumed in [54j which can impact the details of the 
discrimination against the ti background, which is based on a Mt2 analysis. Indeed, 
we find from simulation that the Mt2 variable in our case tends to be rather small, 
and e X A for this analysis is below 0.1% (including the branching ratios). Therefore, 
this search does not set a bound on the RH stop in our scenario. 

There is a second search focusing on fully hadronic top decays [55], that can be seen 
to apply for t^t^ production with — )■ tX^"*" followed by Xi~^ — )■ v^Zjh. For instance, 
when both Z gauge bosons decay invisibly the topology becomes identical to the one 
considered in the above analysis (where — t + ^t)- Also when both Z's decay 
hadronically one has a jet + final state. In fact, although the analysis attempts to 
reconstruct both tops, the required 3-jet invariant mass window is fairly broad. We find 
from simulation that when BR(X|'^ — )■ Zv^ = 1, the ex A of our signal is very similar 
to that in the simplified model considered in [55]. However, when BR(X{'"'" — )■ hug) = 1 
we find that ex A is significantly smaller. Due to the sneutrino vev dependence of these 
branching fractions in our model, we find (for benchmark 1 ) that this search can exclude 
mi^ in a narrow window around 400 GeV for a large sneutrino vev (ve ~ 50 GeV). 
For lower stop masses the search is limited by phase space in the decay — t- tX°"*", 
while at larger masses the sensitivity is limited by the available BR(X|''^ — )• Zue) (see 
Fig. [3|. At small sneutrino vev no bound on m^^ can be derived from this search 
due to the suppressed branching fraction of the Z-channel. We also find that the RH 
stop mass can be excluded in a narrow window around 380 GeV from the decay chain 
Ir — )■ X^~bji followed by — )■ W^u^. Although there are no tops in this topology, 
it is possible for the 3-jet invariant mass requirement to be satisfied, and therefore a 
bound can be set in certain regions of parameter space. 

There is a third search that allows for one hadronic and one leptonic top decay |56|. 
We find that it is sensitive to the LH stop in a narrow window around m^^ ~ 380 GeV 
(for benchmark 1). However, we are not able to set a bound on m^^ from this search. 

We conclude that the present dedicated searches for top squarks are somewhat inefficient 
in the context of our model, but could be sensitive to certain regions of parameter space. 
The most robust bounds on LH stops arise rather from the lepto-quark searches discussed in 
the previous section. However, since the latter do not constrain the RH stop, it is interesting 
to notice that there exist relatively mild bounds (below the top mass) for t^j, as discussed 
above, and perhaps sensitivity to masses around 400 GeV. 



Direct sbottom searches : Ref. 57 sets a limit on the sbottom mass of about 420 GeV, 
based on bb* pair production followed by 6 — )■ txi and Xi ^~Xi i^oi m^o = 50 GeV, and 
assuming BR's = 1). This is essentially our topology when b^ — > tXf and Xf ^ — > W~Ue. 
When kinematically open, these channels indeed have BR close to one, so that the previous 
mass bound would approximately apply (the masslessness of the neutrino should not make an 
important difference). However, BR(6l — > tX^^) can be suppressed near threshold, as seen 
in Fig. 8 For instance, if BR{bL tX^^) = 0.5, the mass bound becomes m^^ > 340 GeV. 
The RH sbottom in our model does not have a normal chargino channel (but rather a decay 
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Signature 


a X e X A[{h] 


L J 


Reference 


> 16 + > 4 jets + 1 lepton + 


8.5 - 22.2 


2.05 


ATLAS 


59 




> 26 + jets + Xt 


4.3-61 


2.05 


ATLAS 


59 




> 36 + jets + Xt 


1.5-5.1 


4.7 


ATLAS 


60 


1 



Table 2: Generic searches for events with b tagged jets. 



1000 



800 



Neutralino "LSP" Scenario 



involving an electron or tau, which falls in the lepto-quark category), so this study does not 
directly constrain fn^^. 

CMS has recently updated their a^- 
based search for sbottom pair produc- 
tion decaying via 6 — j- 6 + Xr [ssj . 
For mLSP = and BR(6 b +Xt) = 
1, they set an impressive bound of 
~ 550 GeV. Taking into account 
the branching fraction for the — )■ 
+ T^T)bR decay mode in our model, 
we find a lower bound that ranges from 
mf''' ^ 330 GeV to m?'^ ^ 490 GeV 
as /i (~ f^xo) ranges from 100 GeV 
to 300 GeV (for our benchmark values 
of the other model parameters). The 
corresponding lower bound on the RH 

470 GeV, in- 



200 



sbottom mass is mp™ 

bR 

dependent of /i. Here we have assumed 
that A333 saturates the bound in Eq. Q, 
as we have been doing throughout. If 
this coupling is instead negligible, thus 
closing the channel, we find that 

490 GeV 



P & 2"'^ Generation 




bL (LQ) 



1 11 (model) 



bL (SUSY search) 



100 



150 



200 
ju [GeV] 



250 



300 



rrir 



290 GeV to mf'^ 

bL 



as /i ranges from 100 GeV to 300 GeV, 
while mP™ ?a 430 GeV, again indepen- 
dent of /i. We note that the bound on 
m^^ sets indirectly, within our model, 
a bound on the LH stop, since the lat- 
ter is always heavier than 6 l (recall that 
the LR mixing is negligible due to the 
approximate i?-symmetry). Typically, 



30 - 50 GeV. 



Figure 14: Summary of exclusions for the "neutralino LSP 
scenario" with Mf = 1 TeV, Mf = 1.5 TcV, Af = and 

= 1, as a function of fi (approximately the X^^ /X^ 
mass). The exclusion on first and second generation squarks 
comes from jets + Xr searches. The bound on bi, come from 
direct bbx'^x'^ SUSY searches, which are somewhat stronger 
than the corresponding lepto-quark searches [dashed line 
marked 5i(LQ)]. This implies, indirectly, a bound on 
about 30-50 GeV larger. (We do not show the bound of 
m^^ ~ 470 GeV, which is independent of /x.) We also indi- 
cate by dashed lines the II lepto-quark searches in the most 
constraining cases: small tan /3 {rrbb search) and large tan/3 
{eejj search). However, these can be completely evaded for 
other values of tan/3. 



Generic searches sensitive to third generation squarks: In Table [2| we summarize a num- 
ber of generic searches with b-tagging and with or without leptons. We see that the bounds 
on a X e X A range from a few fb to several tens of fb. We find that our model cross section 
for these signatures (from pair production of 400 GeV t^, iji, bi or 6^;) are in the same 
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ballpark, although without taking into account efficiencies and acceptance. Thus, we regard 
these searches as potentially very interesting, but we defer a more detailed study of their 
reach in our framework to the future. 



We summarize the above results in Fig. [T4| where we also show the bounds on the first 
two generation squarks (Section|4]), as well as the lepto-quark bounds discussed in Section 5.1 



(shown as dashed lines). The blue region, labeled "6l(SUSY search)", refers to the search 
via two 6-tagged jets plus which has more power than the LQ search that focuses on the 
same final state. The region labeled "^^(model)" refers to the bound on inferred from the 
SUSY search on bi. We do not show the less sensitive searches, nor the bound on 6^, which 
is independent of fi, and about 470 GeV in our benchmark scenario. 



6 Summary and Conclusions 

We end by summarizing our results, and emphasizing the most important features of the 
framework. We also discuss the variety of signals that can be present in our model. Although 
some of the individual signatures may arise in other scenarios, taken as a whole, one may 
regard these as a test of the leptonic i?-symmetry. The model we have studied departs 
from "bread and butter" SUSY scenarios (based on the MSSM) in several respects, thereby 
illustrating that most of the superpartners could very well lie below the 1 TeV threshold in 
spite of the current "common lore" that the squark masses have been pushed above it. 

There are two main theoretical aspects to the scenario: a) the presence of an approximate 
U{1)r symmetry at the TeV scale, and b) the identification of lepton number as the R- 
symmetry (which implies a "non-standard" extension of lepton number to the new physics 
sector). The first item implies, in particular, that all BSM fermions are Dirac particles. A 
remarkable phenomenological consequence is manifested, via the Dirac nature of gluinos, as 
an important suppression of the total production cross section of the strongly interacting 
BSM particles (when the gluino is somewhat heavy). This was already pointed out in Ref. ^5] 
in the context of a simplified model analysis. We have seen here that the main conclusion 
remains valid when specific model branching fractions are included, and even when the gluino 
is not super-heavy (we have taken as benchmark a gluino mass of 2 TeV). We find that: 

• The bounds on the first two generation squarks (assumed degenerate) can be as low 
as 500 — 700 GeV, depending on whether a slepton (e.g. f^) is lighter than the lightest 
neutralino [X^'^ in our notation; see comments after Eq. (js])]. There are two important 
ingredients to this conclusion. The first one is the above-mentioned suppression of 
the strong production cross section. Equally important, however, is the fact that the 
efficiencies of the current analyses deteriorate significantly for lower squark masses. For 
example, the requirements on missing energy and rries (a measure of the overall energy 
involved in the event) were tightened in the most recent jets + analyses (~ 5 fb~ ) 
compared to those of earlier analyses with < 1 fb^^. As a result, signal efficiencies of 
order one (for 1.4 TeV squarks and 2 TeV gluinos in the MSSM) can easily get diluted 
to a few percent (as we have found in the analysis of our model with 700 GeV squarks 
and 2 TeV gluinos). This illustrates that the desire to probe the largest squark mass 
scales can be unduly influenced by our prejudices regarding the expected production 
cross sections. We would encourage the experimental collaborations to not overlook 
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the possibility that lighter new physics in experimentally accessible channels might be 
present with reduced production cross sections. Models with Dirac gluinos could offer 
a convenient SUSY benchmark for optimization of the experimental analyses. It may 
be that a dedicated analysis would strengthen the bounds we have found, or perhaps 
result in interesting surprises. 

It is important to keep in mind that the previous phenomenological conclusions rely 
mainly on the Dirac nature of gluinos, which may be present to sufficient approximation even 
if the other gauginos are not Dirac, or if the model does not enjoy a full symmetry. 
Nevertheless, the presence of the U{1)r symmetry has further consequences of phenomeno- 
logical interest, e.g. a significant softening of the bounds from flavor physics or EDM's 6] 



(the latter of which could have important consequences for electroweak baryogenesis 12 ,61 ). 
In addition, the specific realization emphasized here, where the /^-symmetry coincides with 
lepton number in the SM sector, has the very interesting consequence that: 

A sizable sneutrino vev, of order tens of GeV, is easily consistent with neutrino mass 
constraints (as argued in [l6], [l]; see also Ref. 13 for a detailed study of the neutrino 



sector). The point is simply that lepton number violation is tied to U{1)r violation, 
whose order parameter can be identified with the gravitino mass. When the gravitino is 
light, neutrino Majorana masses can be naturally suppressed (if there are RH neutrinos, 
the associated Dirac neutrino masses can be naturally suppressed by small Yukawa 
couplings). We have also seen that there are interesting consequences for the collider 
phenomenology. Indeed, the specifics of our LHC signatures are closely tied to the non- 
vanishing sneutrino vev (in particular the neutralino decays: Xi'^ — ?■ Zv'e/hue/W' e^, 
or the chargino decay: X^" — )■ W^Ve)- 

This should be contrasted against possible sneutrino vevs in other scenarios, such as 
those involving bilinear i?-parity violation, which are subject to stringent constraints from 
the neutrino sector. Note also that the prompt nature of the above-mentioned decays may 
discriminate against scenarios with similar decay modes arising from a very small sneutrino 
vev (thus being consistent with neutrino mass bounds in the absence of a leptonic U{1)r 
symmetry). In addition, the decays involving a W gauge boson would indicate that the 
sneutrino acquiring the vev is LH, as opposed to a possible vev of a RH sneutrino (see 



e.g. 62 for such a possibility). 

A further remarkable feature -explained in more detail in the companion paper [Ij- is that 
in the presence of lepto-quark signals, the connection to neutrino physics can be an important 
ingredient in making the argument that an approximate U{1)r symmetry is indeed present 
at the TeV scale. In short: 

• If lepto-quark signals were to be seen at the LHC (these arise from the LQD^ "RPV" 
operator), it would be natural to associate them to third generation squarks (within a 
SUSY interpretation, and given the expected masses from naturalness considerations). 
In such a case, one may use this as an indication that some of the A'j33 couplings are 
not extremely suppressed. The neutrino mass scale then implies a suppression of LR 
mixing in the LQ sector. From here, RG arguments allow us to conclude that the three 
Majorana masses, several A-terms and the /i-term linking the Higgs doublets that 
give mass to the up- and down-type fermions (see footnote |2]) are similarly suppressed 
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relative to the overall scale of superpartners given by Msusy, which is the hallmark of a 
U{1) R symmetry. Therefore, the connection to neutrino masses via a LQ signal provides 
strong support for an approximate U{l)fi symmetry in the full TeV scale Lagrangian, 
and that this symmetry is tied to lepton number, which goes far beyond the Dirac 
nature of gluinos. In particular, it also implies a Dirac structure in the fermionic 
electroweak sector, which would be hard to test directly in many cases. Indeed, in 
the benchmark we consider, the lightest electroweak fermion states are Higgsino-like 
and hence have a Dirac nature anyway, while the gaugino like states are rather heavy 
and hence difficult to access. What we have shown is that the connection to neutrino 
masses can provide a powerful probe of the Dirac structure even in such a case. 

The (approximately) conserved i?-charge, together with electric charge conservation can 
impose interesting selection rules (e.g. allowing 2-body decays of the LH squarks, including 
ti, into a state involving an electron but not involving the next lightest chargino, X^'). Of 
course, eventually the approximate i?- symmetry should become evident in the decay patterns 
of the BSM physics. The above lepto-quark signals, and perhaps signals from resonant single 
slepton production fl] that may be present in more general RPV scenarios, can be amongst 
the first new physics signals discovered at the LHC Although by themselves, these may 
admit interpretations outside the present framework, the "L = R" model has a variety of 
signals that provide additional handles. Some of them are summarized below. 

The presence of fully visible decay modes, in addition to those involving neutrinos, may 
give an important handle in the reconstruction of SUSY events. An example is displayed 



in the left diagram of Fig. 15, where one of the squarks decays via q — )■ jX^~ followed by 
X°~ — )■ ejW^ (with a hadronic W), while the other squark gives off missing energy in the 
form of neutrino(s), which can help in increasing the signal to background ratio. Although 
the combinatorics might be challenging, there are in principle sufficient kinematic constraints 
to fully reconstruct the event. 

Perhaps more striking would be the observation of the lepto-quark decay mode of the 



RH strange squark, as discussed at the end of Section 5.1 The pure LQ event (eejj) would 



allow a clean measurement of rrigj^, which could then be used in the full reconstruction of 



"mixed" events involving missing energy, such as displayed in the right diagram of Fig. [T5| 
Furthermore, if the gluinos are not too heavy, associated production of different flavor squarks 
(one being sr) through gluino t-channel exchange, may allow an interesting measurement of 
the second squark mass. Both of these would offer discriminatory power between scenarios 
with relatively light squarks (~ 700 GeV, as allowed by the /^-symmetry) versus scenarios 
with heavier squarks (e.g. > 1 TeV with ultra-heavy gluinos, as might happen within the 
MSSM), by providing information on the scale associated with a putative excess in, say, the 
jets + channel. 

An important possible feature of the present scenario is the presence of flnal states with 
large third-generation multiplicities. We have seen how in the "stau LSP scenario" squark 
pair production can result in flnal states with multiple r's, often accompanied by one or 
more leptons (e or fj,). Although these may not be the discovery modes due to a reduced 
efficiency compared to more standard squark searches, they remain as an extremely inter- 
esting channel to test the present scenario. Similarly, processes such as qq* — )■ jjX^~^X^~ — > 
UtitJ^t^t^ — i- jjrlr^bfibjitLiL, display all the heavy third generation fermions in the final 
state, and it would be extremely interesting to conduct dedicated searches for this kind of 
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Figure 15: Examples of processes with one fully visible decay chain (thus allowing for mass reconstructions), 
while containing a significant amount of missing energy from the second decay chain (that can help for 
triggering and discrimination against backgrounds). The arrows indicate the flow of L — R number. 



topology. Another extremely interesting 'no missing energy" topology arises in the "neu- 
tralino LSP scenario": qq* — )■ jjX^'^X^' — )■ jje^e~[W^W~ . In the lepto-quark sector, 
signals such as bjib*j^ — )■ e'j^eJ^tLiL/T^T^tLtL have not been looked for experimentally, but 
have been claimed to be feasible in Ref. 
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Needless to say, experimentalists are strongly 
encouraged to test such topologies given the expected importance of the third generation in 
connection to the physics of electroweak symmetry breaking. 

Finally, it is important to note that there may be alternate realizations of an approximate 
U{1)r symmetry at the TeV scale. For example, in the realization in which the /^-symmetry 
is identified with baryon number 63 ,64 , the "LSP" decays predominantly to jets, giving 



rise to events with very little missing energy and hence evading most of the current LHC 
bounds. So, these models may hide the SUSY signals under SM backgrounds. A remarkable 
feature of the realization studied in this paper is that fairly "visible" new physics could still 
be present just were naturalness arguments could have indicated. It is certainly essential to 
test such (and possibly other) realizations if we are to address one of the most important 
questions associated to the weak scale: whether, and to what extent, EWSB is consistent with 
naturalness concepts as understood within the well-tested effective field theory framework. 
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A Simplified Model Analysis 

In this appendix we provide details of the interpretation of a number of ATLAS and CMS 
analysis within the simplified models defined for the "neutralino LSP scenario" in Subsec- 
tion iJl 
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A.l Topology (1) : Zv^ 

The LHC searches relevant for this topofogy are: 

• jets + Xr , 

• Z(ll) + jets +Xt , 

• multilepton (> 3/) + jets + (without Z veto). 



We start with the dilepton + 
jets + channel, basing our discus- 
sion on a CMS analysis with 4.98 fb~^ 
{gg production with g — )■ qqx and x 
Z LSPQ) IMS- Bounds are given for 
X = 1/4,1/2 and 3/4. Identifying niq 
with the "gluino mass" , taking rriLSP = 
0, and adjusting the squark mass un- 
til the experimental upper bound on a 
is matched by our theoretical cross sec- 
tion, we find for x = 1/2: ai{mq ^ 
585 GeV) ^ 0.07 pb and a2(mg ^ 
650 GeV) ^ 0.06 pb. What this means 
is that this topology /analysis gives a 



10^ 



lower bound of rrin 



585 GeV when 



Xf"*" is produced as in our scenario, and 
of rrig ^ 650 GeV in a scenario where all 
the squarks decay into neutralino plus 
jet, followed by the decay X^^ — Zu^ 



Z(ll)+ jets + 




V'^SMS K = 100 GeV) 



with BR 

as explained in Subsection 4.2). 
lighter 



[X 



(Ji and (J2 are computed 
For a 

1/4), the correspond- 



bounds are 



ma 



360 GeV and 



mg 

440 GeV, respectively. All of these can 
be read also from Fig. 16 

As a check, and to evaluate the effect 
of the additional two jets in the topol- 
ogy considered in 40 compared to the 



600 
rriq [GeV] 

Figure 16: Production cross-sections of X^^X^~ 
via squark decays, for 2 TeV gluinos (see Subsec- 
tion 4.2 for the definition of cxi and (T2, where cxi 
is the relevant one in our scenario). For reference, 
we show the MSSM total strong production cross- 
section (squarks and gluinos). The dashed lines are 
the SMS upper limit from the CMS searches for the 
channel Z (U) + iets + , assuming m^^ = 100 GeV 
300 GeV |23l. 



and rriy. 



squark pair-production of our case, we have obtained the ex A from simulation of our signal 
{qq production with q — )• qX^^ and X^^ — )• Zue, taking m-^o = 200 GeV) in the various 
signal regions of the CMS analysis]^ We find that the strongest bound arises from the 
"MET Search" with > 300 GeV (with e x A ^ 1.5%, including the branching fractions 
of the Z), and corresponds to a model cross-section of about 40 fb. This translates into the 



^^Note that this topology is not identical to ours, having two extra jets. 

"'^''We have also simulated the case of gg production with g — >■ qqx and x — > Z LSP, taking 900 GeV 
gluinos, heavy (5 TeV) squarks, a massless LSP and x = 1/2, i.e. = 450 GeV. We reproduce the e x A 
in [40l within 30%. 
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Search 


cr X e X A [fb] 


C [fb ^] 


Reference 


1 lepton 


1.1 - 1.7 

1 — 9 


5.8 

4 7 


ATLAS 
ATLAS 


43 
42 




2 OS leptons 


1-5 


4.98 


CMS 


65 




2 bb leptons 


1.6 


2.05 


ATLAS 


66 




z{i+i-) 


0.6-8 


4.98 


CMS 


40 




Multilepton 


1.5 (no Z), 3.5 {Z) 
1 - 2 


2.06 
4.7 


ATLAS 
ATLAS 


41 
42 





Table 3: Upper limits on a x e x A for a number of leptonic channels, with the corresponding luminosity 
and the ATLAS or CMS reference. 



bounds rriq > 640 GeV (based on ai) and rriq > 690 GeV (based on (72), which are somewhat 
stronger than above. 

For the jets + signal we use an ATLAS search with 5.8 fb^ which includes five 

different signals regions depending on the jet multiplicity. In order to apply this analysis, 
we estimate the efficiency times acceptance in our model in the different signal regions 
by simulating our signal (X°"*'X{'~jj production via the processes defining ai, followed by 
X^"*" — )■ with BR = 1), and then applying the cuts in Il9]. Our topology, and our model in 
general, is distinguished by long cascade decays, and we find that the strongest bound arises 
from signal region D (tight), i.e. a 5 jet region, setting a bound on the signal cross-section 
of about 20 fb. We find a lower limit of rriq ~ 635 GeV (based on cti), and rriq ~ 685 GeV 
(based on (72). We conclude that the bounds from this analysis are very comparable to 
those from the Z{ll) + jets + channel. We note that CMS has a MT2-based Simplified 



Model analysis of the jets + signature with 4.73 fb ^ 21 (SMS: gg production with 
g —> qq + LSP). Applying the procedure detailed at the end of Subsection 4.2 we find that 
aiiniq ^ 350 GeV)xBR(Z ^ jj)^ ^ 0.9 pb and a2(mq ^ 440 GeV) xBR(Z jj^ ^ 0.4 pb. 
The fact that these limits are much weaker than those obtained from the ATLAS study may 
be related in part to the additional hard jets that differentiate the gluino from the squark 
pair production topology. 

There are no SMS limits on multilepton searches applicable to our topologies, but there 
are a number of model-independent upper bounds on cr x e x A, as summarized in Table [3] 
However, putting in the BR(Z — t- and taking into account the general lessons from the 
computed efficiencies for "Topologies (3) and (4)" below, we conclude that such searches are 
less sensitive than the previous two searches. 



A.2 Topology (2) : Xf+ hue 

For this topology we use the ATLAS jets + search since the Higgs decays mostly 
into hadrons. This topology is characterized by a high jet multiplicity, as was the case with 
the hadronic Z of the previous topology. The efficiency times acceptance is the same as in 
the case studied above (with a Z instead of h), so that the bound on the model cross section 
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is about 20 fb. We find a lower limit of rriq ~ 605 GeV (based on ai), and rriq ~ 655 GeV 
(based on (72). 

Note that, since the Higgs decays predominantly into bb, searches with b tagged jets 
are interesting for this topology. In a search for final states with and at least three 
b-jets (and no leptons), ATLAS sets a bound on the corresponding visible cross section of 
about 2 fb |60]. However, simulating our signal (for 700 GeV squarks) in MG5 + Pythia + 
Delphes, we find an extremely small efficiency for the present topology: e x A ^ 10~^ for 
their signal regions SR4-L and SR4-M (and much smaller efficiencies for the other SR's). 
This arises from the aggressiveness of the requirement and the combined efficiency of 
tagging three 6-jets. As a result we infer a very mild bound on the model cross section of 
about 18 pb, and no meaningful bound on the squark masses, as such a cross section can be 
reached only for squarks as light as a couple hundred GeV, where the e x A would be even 
smaller. Nevertheless, it would be interesting to optimize such an analysis for the present 
model (with suppressed production cross sections), and furthermore to try to reconstruct bb 
resonances at about 125 GeV. 

The leptonic searches are not constraining due to the significant suppression from the 
Higgs branching fraction into final states that might involve leptons. 



A.3 Topology (3) : W-el 

In this case the two relevant searches are: jets + two leptons without ^t, and multileptons 
+ jets + The first signal has a branching fraction of BR(W^ jjY ~ 0.45. However, at 
the moment there are no searches that constrain this topology, since these typically include 
important cuts on the missing transverse energy. It would be interesting to perform a 
dedicated search for this signal. Here we focus on the existing multilepton searches. CMS 



has a detailed analysis including a large number of channels 67 . Unfortunately, the results 



are model-dependent and no information on o" x e x A for the different signal regions is 
provided. ATLAS has a > 4 leptons ( + jets + ^t) search with and without Z veto 41 



Their upper limit (with a Z veto) is a x e x A ^ 1.5 fb. We find from simulation of our 
signal that, for this analysis, ex Ak. 0.02 (which includes the branching fractions of the W 
decays). We can therefore set a limit of rriq > 580 GeV (based on cxi) and rriq > 630 GeV 
(based on (72), corresponding to a model cross section of about 75 fb. 



A.4 Topology (4) : X+ W+Ue 
In this case the relevant LHC searches are: 

• jets + Xr , 

• 1 lepton + jets + , 

• OS dileptons + jets + ■ 

^^This corresponds to combining a number of channels with different flavor composition, not all of which 
are present in our model. Thus, this result provides only an estimate for the possible bound in our model 
from such a multi-lepton search. 
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We can use again the ATLAS bound on jets + discussed above. In this case, however, 



the efficiency times acceptance turns out to be smaller, The strongest constraint arises 
again from signal region D (tight) in [l9j, and gives an upper bound on our model cross 
section of about 60 fb. This translates into rrig > 530 GeV (based on cxi) and nig > 650 GeV 
(based on (T2). 

In a multi-lepton study, the ATLAS collaboration has considered our simplified model 
(model C in [42]), except that all the squarks are assumed to decay with unit branching 
fraction through the chargino channel (i.e. the process characterized by (J2). If we assume 
the same efficiency times acceptance for the process in our scenario, i.e. based on ai, we read 
from their Fig. 10 the bounds rriq > 410 GeV (based on ai), and niq > 500 GeV (based on 
(T2), which correspond to model cross sections of about 0.35 pb. 

ATLAS also has a search for exactly 1 lepton + > 4 jets + ^t, setting a bound on 
ax ex A ^ 1.1 — 1.7fb, depending on whether the lepton is an electron or a muon f43|. 
Simulation of the above process {qq production with q — )• qX^~ followed by — )• W^Vei 
taking m-^± = 200 GeV) gives that the ATLAS analysis has ex Ak, 10~^ (this includes the 
branching fractions for the W decays). We see that the efficiency is quite low. This is due, 
in part, to the fact that the analysis requires at least four jets with pt > 80 GeV. While the 
two jets from squark decays easily pass the pt cut, the other two jets arise from a W decay 
and are softer (the other W decaying leptonically). But when the quarks are sufficiently 
boosted to pass the pr cut, they also tend to be coUimated, and are likely to be merged into 
a single jet. As a result, using an upper bound on the model cross section of order 1 pb, we 
get a rather mild bound of rriq > 350 GeV (based on ai), and nig > 430 GeV (based on 0-2) • 
For the OS dilepton signal CMS sets a bound ofcrxexA< l-5fb with 4.98 fb~^ [65]. 
Our simulation gives e x A ~ 10~^ (including the W BR's), resulting again in an upper 
bound on the model cross section of about 1 pb, and the same mild bounds as above. 

CMS studies a simplified model {gg production with gi — )• 92 — > (IIX^) with 



4.98 fb ^ 23 , where the neutralino is the LSP, while decays into W^x^- Therefore, as 



in our scenario, a single lepton is produced via W decay, although there are two extra hard 



jets compared to our case from the gluino versus squark production. From Fig. 8 of 23 , with 
?^LSP = 0, we find that our cross section, in the range 300 GeV < rriq < 800 GeV, is more 
than an order of magnitude below the current sensitivity. Here we used our cxi including the 
branching for exactly one of the W's to decay leptonically. 

Finally, our model has very suppressed SS dilepton signals due to the Dirac nature of 
the gluinoj^ so that no interesting bounds arise from this search. In conclusion, for this 
simplified topology, the strongest bounds again arise from the generic jets + searches, 
although it should be possible to optimize the leptonic searches to our signal topologies to 
obtain additional interesting bounds. 



^^From simulation via MG5 + Pythia + Delphes of ^jj via the processes in the definition of ai 



(see Subsection 4.2). 

^'^Two SS positrons can be obtained tlirough ulUr production consistent with the Dirac nature of gluinos, 
but the SS dilepton cross section is small, at the 0.2 fb level for 700 GeV squarks in the "neutralino LSP" 
scenario. In the "stau LSP" scenario the SS dilepton + jets + signal can reach 1 — 2 fb. 
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